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ABSTRACT 

The  differential  equations  and  boundary  conditions  describing  the 
periodic-flow  heat  exchanger  (frequently  referred  to  as  rotary  regen- 
erator) are  sufficiently  complicated  to  preclude  complete  analytical 
solution.  Because  of  the  interest  in  this  type  of  heat  exchanger  con- 
siderable effort  has  been  expended  to  obtain  useful  solutions  to  the 
analytical  problem  of  predicting  its  performance.  Several  approximate 
solutions  have  been  proposed  but  each  of  these  Is  limited  due  to  the 
restricted  ranges  of  parameters  in  which  they  apply. 

A  numerical  finite-difference  method  of  calculating  the  spatial 
temperature  distribution  in  the  periodic -flow  type  heat  exchanger  is 
presented  by  considering  the  metal  "stream"  in  cross -flow  with  each  of 
the  fluid  streams  as  two  separate  but  dependent  exchangers.  In  the 
development  no  assumptions  are  made  which  would  restrict  the  range  of 
parameters  in  which  the  analysis  would  be  applicable.  The  exchanger 
effectiveness  has  been  evaluated  over  the  following  ranges  of  dimen- 
sionless  parameters: 

i*0   ^  Cmln/Cmax  ^  0.10 

1.0  ^  Cr/Cmin  ^  °° 

1.0  ^   NTU0  ^   10 

1.0  ^   (hA)«  ^   0.25 
The  National  Cash  Register  102A  general  purpose  digital  computer 
was  used  to  accommodate  the  large  number  of  subdivisions  necessary  for 
accuracy  and  an  extrapolation  of  data  to  zero  element  area  (infinite 
subdivision)  was  used  to  obtain  values  of  effectiveness  good  to  four 
significant  figures. 

The   author  expresses  his  appreciation  to  Assistant  Professor  C.P. 
Howard  for  his  direction  and  encouragement  in  this  work. 
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NOMENCLATURE 
SYMBOLS 

A  =  heat  transfer  area  on  side  designated  by  subscript,  sq  ft 

C  =  heat  capacity  rate  (Wc)  of  fluids  or  rotor  matrix  according  to 

subscript,  Btu/hr,  deg  F 
c  =  specific  heat  of  fluids  (constant  pressure)  or  rotor  matrfx 

material  depending  on  subscript,  Btu/lb,  deg  F 
h  =  unit  conductance  for  thermal  convection  heat  transfer, 

Btu/hr,  sq  ft,  deg  F 
k  =  unit  thermal  conductivity,  Btu/hr,  sq  ft,  deg  F/ft 
K  =  a  constant,  defined  as  used  and  distinguished  by  numerical 

subscript 
N  =  number  of  subdivisions  of  fluid  or  matrix  solid  flow  "stream" 

according  to  subscript 
Q  =  heat  transfer  rate,  Btu/hr 
T  =  temperature  of  fluid  or  matrix,  inlet  or  outlet,  depending  on 

subscripts,  deg  F 
W  =  mass  flow  rate  of  fluid  (lb/hr)  or  matrix  (rev/hr)(mass) 

according  to  subscript 
DIMENSIONLBSS  PARAMETERS 

E  -  exchanger  heat  transfer  effectiveness;  ratio  of  actual  to 

thermodynamically- limited  maximum  possible  heat  transfer  rate 
Cmin/Cmax  =  capacity  rate  ratio  of  fluid  flow  streams 
Cr/Cmin  =  capacity  rate  ratio  of  rotor  matrix  to  minimum  fluid 

capacity  rate 
(hA)'  =  (hA)n/(hA)x,  geometric  symmetry  and  heat  transfer  balance 

of  exchanger 


NTU  =  (hA)/C,  number  of  transfer  units  on  side  designated  by 
subscript 


«TOo=   -7T—  /  //„.,  .  / 


over-all  number  of  transfer  units 
SUBSCRIPTS 

avg  «  average 

r   =  rotor  matrix 

i   =  inlet 

o   =  outlet,  or  over-all  in  case  of  NTU0 

min  *  minimum  magnitude 

max  =  maximum  magnitude 

n   =  value  on  side  of  C^q 

x   =  value  on  side  of  C^qx 

c   =  cold  side 

h   =  hot  side 

oo  =  subscript  on  E  to  indicate  value  when  CT/Cmin  =  oO 
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1.  Introduction 

In  the  past  decade  considerable  interest  has  been  shovn  in  the  ap- 
plication of  the  periodic-flow  heat  exchanger  (rotary  regenerator)  to  the 
Gas-turbine  cycle  for  exhaust-gas  thermal- energy  regeneration.  Regen- 
eration by  whatever  method  applied  to  the  gas-turbine  cycle  obtains,  among 
other  improvements,  higher  thermal  efficiency,  lover  specific  fuel  con- 
sumption, Improved  part  load  performance  and  lower  optimum  pressure  ratio. 
Particular  advantages  to  be  expected  from  the  application  of  the  periodic- 
flow  heat  exchanger  to  the  gas-turbine  process  are  discussed  by  Harper 
and  Rohsenow  (l)1  and  Coppage  and  London  (2),  among  others.  Principal 
among  these  is  the  possibility  of  a  more  compact  unit  than  is  possible 
with  the  conventional  stationary  surface  types.  This  is  due  to  the  fact 
that  the  periodic-flow  heat  exchanger  matrix  can  be  finely  divided,  with- 
out structural  difficulty,  resulting  in  an  appreciably  reduced  heat  ex- 
changer volume  relative  to  a  stationary  type  for  a  given  effectiveness 
and  pressure  drop.  The  major  disadvantage  of  the  rotating  periodic-flow 
type  is  the  sealing  difficulties  encountered  and  the  attendant  leakage 
involved  when  an  appreciable  pressure  ratio  is  used.  This  has  been  re- 
duced to  acceptable  limits  in  at  least  two  current  development  programs 
(3,U)  and  further  refinements  should  follow. 

Because  of  the  interest  in  this  type  of  heat  exchanger  there  has 
been  considerable  effort  expended  to  obtain  useful  solutions  to  the 
analytical  problem  of  predicting  its  performance.  The  differential 
equations  and  boundary  conditions  describing  the  system  are  sufficiently 
complicated  to  preclude  complete  analytical  solution.  An  analytical 


•'■Numbers  in  parentheses  refer  to  the  Bibliography  at  the  end  of  the  thesis. 


solution  is  available  for  the  special  case  of  infinite  matrix  rotative 
speed  where  the  expression  for  effectiveness  takes  the  same  form  as  that 


of  a  counter-flow  direct  type  (2)  and  is  given  by 

r  ._       /  -  e r  7 


and  for  the  further  condition  of  C^jj/Cj^^  ■  unity,  reduces  to 

p  =     w*  [uJ 

Consequently,  several  approximate  methods  of  solution  have  been 
proposed.  Perhaps  the  oldest  is  that  of  Hausen  (6),  described  and  elab- 
orated by  Johnson  (8,9),  which  begins  with  the  known  solution  for  the 
initial  cooling  of  a  matrix  element,  initially  at  uniform  temperature, 
by  the  passage  of  a  cold  fluid.  By  dividing  the  regenerator  tube  into  a 
large  number  of  strips  of  equal  length  and  applying  this  elemental  solu- 
tion to  each  strip  in  turn,  a  family  of  simultaneous  linear  algebraic 
equations  may  be  derived  whose  solution  gives  the  desired  temperature 
variation  in  the  periodic  case. 

A  second  method,  due  to  Iliffe  (10),  uses  the  solution  to  the  prob- 
lem of  the  initial  cooling  of  a  matrix  element  at  arbitrary  initial  temp- 
erature distribution,  derived  by  Nusselt.  Using  this  elemental  solution, 
the  method  of  Iliffe  leads  to  two  simultaneous  integral  equations  which 
require  numerical  methods  for  their  solution.  In  a  third  method,  Saun- 
ders and  Smoleniec  (7)  replace  the  governing  partial  differential  equations 
by  their  equivalent  finite  difference  equations  and  effect  an  approximate 
solution  by  relaxational  techniques.  In  a  fourth  approach,  Tipler  (l6) 
replaces  the  thermal  problem  by  its  electrical  analogue  and  observes  the 
solution  on  an  oscilloscope.  In  yet  a  fifth  approach,  the  perturbation 


method  used  by  Schultz  (8)  has  been  extended  by  Jones  and  Fax  (5)  lead- 
ing to  algebraic  statements  for  the  departure  of  the  regenerator  effec- 
tiveness from  its  asymptotic  value,  given  by  equation  [ij   ,  for  several 
types  of  asymmetry.  Each  of  these  solutions  is  limited  in  one  respect 
or  another,  primarily  due  to  the  restricted  ranges  of  parameters  in  which 
they  apply. 

It  is  the  purpose  of  this  thesis  to  present  the  results  of  a  finite- 
difference  numerical  analysis  which  was  carried  out  with  the  aid  of  a 
digitial  computer  such  that  there  were  no  limitations  imposed  on  the  four 
independent  dimensionless  parameters  selected. 
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2.  Method 

The  conventional  Idealizations  and  boundary  conditions  assumed  in 
the  usual  derivation  of  the  governing  differential  equations  (2)  are  as 
follows : 

a.  The  thermal  conductivity  of  the  matrix  is  zero  in  the  direction 
of  fluid  flow  and  matrix  metal  "flow".  It  is  infinite  In  the 
other  direction  normal  to  the  fluid  flow. 

b.  The  specific  heats  of  the  two  fluids  and  matrix  material  are 
constant  with  temperature. 

c.  No  leakage  of  the  fluids  occurs  either  due  to  direct  leakage 
or  carry  over,  and  each  fluid  flow  is  unmixed. 

d.  The  convective  conductances  between  the  fluids  and  the  matrix 
are  constant  with  flow  length. 

e.  The  fluids  pass  In  counterflow  directions. 

f .  Entering-fluid  temperatures  are  uniform  over  the  flow  inlet 
cross  section  and  constant  with  time. 

g.  Regular  periodic  conditions  are  established  for  all  matrix 
elements,  i.e.,  steady  state  condition. 

With  these  idealizations  a  typical  element,  Fig.  1,  from  each  fluid 
side  of  the  exchanger  can  be  represented  schematically  as  shown  in  Fig. 
2.  Although  Fig.  1  is  a  disk  (axial-flow)  type  exchanger,  typical  el- 
ements of  a  drum  (radial- flow)  type  can  be  selected,  with  the  same  ass- 
umptions, which  can  be  represented  schematically  in  the  same  manner  as 
Fig.  2. 

It  Is  apparent  from  Fig.  2  that  each  of  these  elements  can  be  regard- 
ed as  a  crossflow  exchanger  with  a  gas  stream  and  a  metal  "stream".  This 
extension  of  the  method  of  calculating  crossflow  exchangers  to  the  rotary 
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regenerator  vas  suggested  by  G.  M.  Dusinberre  in  his  discussion  of  ref- 
erence (2).  It  should  be  noted  that  the  Trl  and  Tr0  on  the  two  elements 
are  not  the  same,  but  to  avoid  further  subscript  notation  the  side  in 
question  will  be  specified  when  it  is  not  obvious.  However,  the  Tro  of 
a  particular  element  is  the  Tr^  of  the  next  adjacent  element  In  the 
direction  of  matrix  flow,  and  this  is  of  particular  significance  at  the 
seals  where  the  direction  of  fluid  flow  reverses. 

Consider  for  the  moment  only  the  heat  exchange  element  on  the  side 
of  Cjuax,  Fig  2(a).  For  simple  unmixed  cross  flow  with  uniform  tempera- 
tures across  the  inlets  and  average  temperatures  at  the  outlets,  the 
heat  transfer  rate  for  this  element  is 

Q  =  <WTXi  -  Txo)(l/Nx)     [2] 

Q  =  Cr(Tro  -  Trl)(l/Nr)       [3] 

Q  =  (hA)x  ATavg(l/NrNx)     M 

For  a  small  enough  element  the  arithmetic-mean  temperature  differ- 
ence may  be  assumed  valid,  so  that 

ATavg  =  (l/2)(Txl  +  Txo)  -  (l/2)(Trl  +  Tro)   .  .  .  .  [5] 
By  using  equations  [2]   ,  [3]  ,  [k]   ,  and  /5J  the  outlet  temperatures m 
of  the  two  streams  may  be  solved  for  in  terms  of  the  inlet  temperatures 
and  put  in  the  form  (Appendix  1) 

Txo  -  Txl  -  K,(Txl  -  Trl)    [6a] 

Tro  "  Tri  +  lfe(Txl  -  Tri)    [6bJ 

By  the  same  method  the  outlet  temperatures  of  the  element  on  the 
side  of  Cain  may  be  solved  for  in  terms  of  its  inlet  temperatures  giving 

Tno  -  Tni  +  K3(Trl  -  Tni)    [?aj 

Tro  =  Tri  -  K4(Trl  -  Tni)    [7b] 

By  using  the  definitions  of  the  parameters  and  proper  algebraic 


manipulation  the  foregoing  constants  may  be  put  directly  in  terms  of 

the  dimensionless  parameters  and  the  number  of  subdivisions  of  the  matrix 

and  fluid  streams,  giving 
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Now  consider  the  schematic  representation  of  the  composite  heat 
exchanger  made  up  of  elements  as  shown  in  Fig.  5.  For  illustration  the 
fluid  streams  and  the  matrix  stream  have  each  been  divided  into  three 
equal  substreams  to  form  the  elements.  The  double  lines  represent  the 
area  of  the  seals.  It  should  be  noted  that  the  left  edge  is  physically 
the  same  as  the  right  edge  and  therefore  the  matrix  inlet  temperature  of 
a  particular  matrix  substream  on  the  left  must  be  identical  to  the  matrix 
outlet  temperature  of  that  substream  on  the  right.  This  will  be  referred 
to  later  as  the  reversal  condition. 

In  Fig.  3  where  individual  element  temperatures  are  designated  with 
the  same  label  it  is  not  implied  that  they  are  equal  numerically,  except 
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where  indicated  as  constant  across  the  fluid  Inlets.  This  is  done  to 
avoid  additional  subscript  notation.  Also,  it  is  not  to  be  implied  from 
Fig.  3  that  there  is  necessarily  any  symmetry  with  respect  to  the  areas 
involved . 

Now  for  calculation  purposes  if  the  fluid  Inlet  temperatures  are 
given  some  convenient  values,  such  as  one  and  zero,  and  a  matrix  inlet 
temperature  distribution  is  assumed  on  the  left  edge,  the  remaining 
temperatures  of  each  element  can  be  calculated  in  the  order  indicated 
by  the  element  numbers  by  repetitive  use  of  equations  [6]  and  then  jlj   . 
If  the  temperature  distribution  assumed  on  the  left  edge  was  correct 
then  it  would  be  duplicated  on  the  right,  i.e.,  the  reversal  condition 
was  fulfilled,  and  the  problem  would  be  solved  for  the  particular  set 
of  parameters  used.  If,  however,  this  is  not  the  case,  then  the  result- 
ing temperature  distribution  on  the  right  is  now  used  on  the  left  and  the 
procedure  repeated  until  the  reversal  condition  is  met. 

After  several  iterations  of  this  type  the  cyclic  character  of  the 
calculations  for  each  element  and  column  on  a  side  becomes  apparent  and 
the  beauty  of  automatic  computing  equipment  is  fully  appreciated.  It 
can  be  seen  that  If,  in  the  process  of  computing  the  outlet  temperatures 
of  element  one,  the  Tx0  is  put  where  the  Txl  was  and  a  control  number 
increased  by  one  so  that  the  Tr^  of  element  two- will  be  used  next  the 
same  sequence  of  operations  can  now  be  used  to  compute  element  two,  and 
so  on  down  the  first  column  of  elements.  Also,  if  the  Tro  of  each  ele- 
ment is  put  where  its  Trl  was,  then  upon  reaching  the  end  of  the  column 
and  "zeroing"  the  above  control  number  the  same  cyclic  operation  can  be 
used  for  the  next  column.  A  second  control  number  is  used  to  indicate 
when  one  side  is  finished  and  then  a  similar  procedure  is  followed  on 
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the  other  side.  In  this  scheme  it  is  only  necessary  to  provide  stor- 
age for  the  TTi   and  Tj^  of  each  row,  which  will  be  used  for  comparison 
purposes  to  determine  if  the  reversal  condition  is  met,  and  the  Txo  of 
each  column  on  one  side  and  the  Tno  of  each  column  on  the  other  side, 
which  will  be  used  to  make  a  heat  balance  and  compute  the  effective- 
ness when  the  reversal  condition  is  met. 

That  the  reversal  condition  eventually  will  be  met  by  this  iterative 
scheme  can  be  shown  (Appendix  1)  to  depend  on  the  physical  conditions  of 
the  problem  and  the  second  law  of  thermodynamics.  This  will  be  called 
the  condition  for  convergence  and  for  this  problem,  using  equal  sub- 
divisions of  each  stream,  it  may  be  stated  in  the  following  form: 

NTU0["l  +  l/(hA)'J  Cmin/Cmax  ^  2N [aj 

Accepting  for  the  moment  the  idealizations  and  assumptions  made, 
then  the  accuracy  of  the  solution  obtained,  by  the  very  nature  of  a 
finite-difference  process,  depends  only  on  the  number  of  subdivisions 
used.  It  is  apparent  from  condition  \s]   that  a  greater  number  of  sub- 
divisions will  also  enhance  the  convergence. 

A  considerable  number  of  subdivisions  can  be  handled  by  automatic 
computing  equipment  to  obtain  a  certain  accuracy  and  aid  convergence, 
but  the  number  of  subdivisions  used  must  be  arrived  at  by  a  compromise 
between  the  accuracy  desired  and  time  available.  The  easiest  way  to 
illustrate  the  accuracy  obtained  for  a  given  degree  of  subdivision  is 
to  solve  the  problem  for  several  values  of  subdivision  and  extrapolate 
the  results  obtained  to  zero  element  area.  As  an  example,  consider  the 
following  results  for  the  indicated  parameter  values: 

Cmin/Cmax  -  0.90,  NTU0  =  6,  Cr/Cmin  =  (hA)'  =  1.0 
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No.  subdivisions    Effectiveness 
8  0.78l6 

16  0.7800 

32  0.7796 

These  results  are  plotted  versus  element  area  and  extrapolated  to 
zero  In  Fig.  3a  where  it  Is  found  that  the  Indicated  effectiveness  is  0.7795- 
This  correction  is  not  constant  for  all  ranges  of  parameters,  being 
largest  for  low  values  of  Cr/Cmin  and  high  values  of  NTU0  ,  the  fore- 
going example  being  typical.  The  procedure  used  in  this  work  was  to 
calculate  Tables  1  through  9  with  16  subdivisions,  then  recalculate 
Tables  1  and  9  with  32  subdivisions.  Table  9  was  calculated  to  show 
the  effect  of  varying  the  parameter  (hA)1,  and  it  also  contains  rep- 
resentative values  from  the  other  tables  so  that  with  these  data  all 
of  the  tabulated  data  were  extrapolated  to  zero  element  area. 

The  computer  used  has  the  equivalent  of  nine  significant  decimal 
figure  accuracy;  it  was  specified  that  the  reversal  condition  be  ful- 
filled to  only  five  figures.  Of  more  significance  than  this  was  the 
final  requirement  that  before  a  solution  was  considerd  as  being  obtained 
the  heat  balance  error,  which  is  a  direct  measure  of  error  in  the  effec- 
tiveness introduced  by  accuracy  of  computation,  be  less  than  0.024  per 
cent. 

The  remainder  of  the  computer  program  was  essentially  a  system  of 
"feeding"  it  the  next  set  of  parameters  after  each  solution.  Appendix 
2  describes  the  computer  program  and  its  various  modifications  in  some 
detail.  Several  interesting  features  incorporated  in  it  are: 

a.  Up  to  6k   subdivisions  could  be  accommodated  without  external 
memory  storage  simply  by  modifying  certain  control  commands. 
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b.  Final  matrix  temperature  distributions  were  punched  on  cards 
for  future  use. 

c.  A  logical  choice  of  initial  estimate  for  matrix  temperature 
distribution  for  each  nev  problem  was  made  from  data  available 
from  work  up  to  that  point. 

d.  Optional  read  out  of  all  or  selected  temperatures  was  provided, 


Ik 


3.  Comparison  of  Results 

The  results  obtained  are  reported  in  tabular  (Tables  1-9)  and 
graphical  form  (Figs,  k   -  11).  Several  general  observations  can  be  made 
from  the  data: 

a.  As  Cmin/Cmax  decreases,  for  given  values  of  the  other  parameters, 
the  effectiveness  increases.  This  is  predicted  by  the  limiting  solution 
for  Cp/C^jj  =  00   ,  equation  [l]   ,  and  vas  pointed  out  in  the  illustrative 
problem  of  reference  (2).  This  is  illustrated  in  Fig.  12. 

b.  By  using  NTU0  as  defined  in  the  nomenclature,  vhich  contains 
(hA)'  explicitly,  the  influence  of  (hA)'  at  a  given  NTU0  is  very  small 
for  CjQjLn/Cmax  ^  0.9  .  This  conclusion  vas  also  reached  by  Coppage  and 
London  (2)  in  their  evaluation  of  the  work  by  Hausen  (6),  Saunders  and 
Smoleniec  (7),  J.E.  Johnson  (8),  and  Iliffe  (10).  For  Cminfcmax   ^  °*9> 
however,  the  influence  of  (hA)1  becomes  increasingly  pronounced,  and  the 
variation  for  (hA)*  =  unity  and  0.25  may  be  as  much  as  seven  per  cent  for 
Cmin/Cmax  =0.1.  See  Table  9. 

c.  The  difference  in  effectiveness  for  Cr/Cmln  ^   5  and  the  value 
given  by  the  limiting  solution  for  Cr/C^n  =  *  is  very  small.  This 
may  be  seen  from  inspection  of  Tables  1  through  8  or  Figs,  h   through  11. 

R.W.  Johnson  (11)  used  much  the  same  method  employed  here  but  had 
limited  access  to  a  computer  and  was  therefore  able  to  calculate  only  a 
relatively  small  number  of  the  points  reported  here.  Also  due  to  limited 
computer  time,  at  most  ten  and  in  some  cases  only  five  subdivisions  were 
used  with  no  extrapolation  to  zero  element  area  correction.  Allowing  for 
this  the  corresponding  results  are  practically  identical. 

Coppage  and  London  (2)  recommended  the  results  of  J.E.  Johnson  (9) 
which  were  obtained  for  Cmin/Cj^^  =  unity.  These  results,  in  terms  of 
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the  parameters  used  here,  were  later  Incorporated  in  reference  (12)  with 
extrapolation  'to  Cmin/Cmax  =  0.7  .  The  variation  in  corresponding  values 
of  Table  1  of  this  work  and  Table  6,  page  11-32,  of  reference  (12)  which 
are  for  C^q/Cj^^  *  unity  is  less  than  one  per  cent.  Only  when  Cmin/Cmax 
=  0.7  and  Cr/Cmin  =  unity,  where  the  extrapolation  method  used  was  ad- 
mitted to  be  most  uncertain,  is  the  variation  slightly  more  than  two  per 
cent. 

An  expression  for  effectiveness  in  closed  form  may  be  desir~able  for 
certain  uses  such  as  computer  analysis  of  plant  performance.  A  relative- 
ly simple  empirical  formulation  is  suggested  in  reference  (12)  for 
Cmin/Cmax^  0.7  and  Cr/Cmin  ^  2  and  is  of  the  form 

*tr  "■  ^  '-  *• (Cr/Cmin)m  •  •  •  — fl0] 

where  Kg  =  l/9  ,  m  =  -2,  and  E^  is  determined  from  equation  j[lj  .  Such 
a  form  cannot  be  applicable  for  all  values  of  parameters  with  single  val- 
ues of  Kg  and  m,  but  it  agrees  reasonably  well  with  the  results  reported 
here  for  the  ranges  specified.  It  has  been  found  that  with  slight  var- 
iation of  Kg  and  m  very  close  agreement  can  be  obtained  in  a  small  range 
of  parameters  which  may  be  of  particular  interest.  As  an  example,  in 
the  range 

l-°  ^  Cmin/Cmax  ^  °-90 
1.25  ^  VC^  ^   5.0 

3.0  ^   NTUQ   ^    9.0 

with  K5  =  1/9  and  m  =  -I.87,  the  error  is  less  than  one  per  cent.  Fig. 

* 
13  illustrates  the  method  used  to  determine  Kg  and  m  in  this  range. 
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14-.  Limitations 

The  most  questionable  of  the  idealizations  made  is  that  of  zero 
matrix  conductivity,  particularly  in  the  direction  of  fluid  flow.  Al- 
though the  effect  of  matrix  conductivity  could  be  included  in  the  finite- 
difference  equations  the  additional  degrees  of  freedom  would  make  adequate 
coverage  of  all  the  parameters  extremely  voluminous  and  time  consuming. 
If  for  a  particular  exchanger  design  the  matrix  conductivity  was  consider- 
ed important  then  the  problem  could  be  solved  individually  on  a  computer 
with  adequate  capacity.  A  simpler  alternative  is  to  neglect  the  conduct- 
ion and  correct  the  effectiveness  for  this  effect.  Prom  theoretical  con- 
siderations Professor  A.L.  London  has  shown  (13)  that  a  correction  factor 
for  conductivity  in  the  direction  of  fluid  flow  for  approximately  equal 
fluid  capacity  rates,  of  the  form 


4E   _   kA 


M 


E      LC 

should  result  in  a  somewhat  pessimistic  prediction  of  the  reduction  in 
effectiveness.   (A  in  this  case  is  the  matrix  cross- sectional  area  avail- 
able for  heat  conduction  in  the  direction  of  fluid  flow  on  both  sides, 
and  L  is  the  flow  length).  Schultz  (lU)  has  solved  for  the  effectiveness 
as  a  function  of  the  conduction  parameter  in  equation  (llj  for  three  cases. 
In  the  nomenclature  used  here,  these  three  cases  are  for  the  following 
conditions: 

a.     Cmin  =  C,^  ,  Cy/C^n  very  large,  k  finite 

b»     Cmin  =  cmax  »  Cr/Cmin  =  °°      >   k  finlte>  regenerator  sub- 
divided into  sections  in  direction  of  fluid  flow 

c»     cmin  *  cmax  ,   cr/cmin  finite,  k  =  °o     ,  regenerator  sub- 
divided into  sections  in  direction  of  fluid  flow. 
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Prom  the  curves  In  reference  (lU)  vhich  allow  direct  comparison  with 
the  results  presented  here  it  is  found  that  for 

*A   ^   0.02    [12] 

LC 

the  reduction  in  effectiveness  is  at  most  about  one  per  cent.  Addition- 
al work,  not  available  at  the  time  of  writing,  concerning  the  influence 
of  conduction  is  reported  in  reference  (15). 

The  effect  of  leakage  can  be  calculated  separately  resulting  in  a 
correction  to  the  effectiveness  without  leakage.  This  has  been  shown  to 
be  small  for  as  much  ten  per  cent  leakage  (1). 
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5.  Summary  and  Conclusions 

The  numerical  procedure  used  in  this  work  assumes  no  symmetry, 
either  with  respect  to  heat  capacity  rates  of  the  fluids  or  area  ratios 
of  the  two  sides,  and  therefore  is  not  limited  in  this  respect.  The 
ranges  of  parameters  which  have  been  covered  effectively  are 

1«°  ^  cmin/cmax  ^  °-l° 

1.0  £   Cr/Cmin  £  °o 

1.0  ^   NTUq  ^   10 

1.0  ^  (hA)'  ^   0.25 
The  results  are  tabulated  at  frequent  enough  intervals  and  over  suffi- 
cient ranges  of  the  parameters  to  include  the  gas-turbine  regenerator 
application  as  well  as  other  applications,  such  as  intercooling,  where 
the  values  of  the  parameters  might  be  quite  different. 

The  error  introduced  by  the  finite-difference  type  of  solution  has 
been  minimized  by  use  of  sufficient  subdivision  and  further  extrapolation 
to  zero  element  area  so  that  four  significant  figures  in  the  results  are 
justified. 

The  effects  of  deviations  from  the  assumed  idealizations  have  not 
been  included  in  the  results. 

Even  though  the  results  obtained  cannot  be  made  to  fit  a  completely 
general  closed  form,  a  relatively  simple  expression  for  effectiveness 
can  be  obtained  which  will  give  good  agreement  over  a  reasonable  range 
of  parameters. 
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TABLE  1 

PERIODIC 

J-FLOW  HEAT 

1  EXCHANGES 

[  EFFECTIVENESS 

Omln/Omax 

=   1.0 

(hA)« 

=   1.0 

L*      /i-  . .  .»  __      „ 

NTU0 

1 

1.25 

1.5 

2 

3 

5 

1 
00 

0.5 

0.3221 

0.3282 

0.3301+ 

0.3329 

0.3333 

1.0 

0.1*665 

0.1*779 

0.1*81*5 

0.1*912 

0.1*960 

0.1*985 

0.5000 

1.5 

0.5W 

0.5757 

O.586I 

0.5978 

0.6000 

2.0 

0.6007 

0.6231 

O.6360 

0.61*91 

O.6587 

O.6638 

0.6667 

2.5 

O.6385 

0.6792 

0.691*1 

0.7110 

0.711*2 

3.0 

0.6672 

0.6959 

0.7119 

0.7280 

0.71*00 

0.7^63 

0.7500 

3.5 

O.69OO 

0.7377 

0.75^6 

0.7739 

0.7778 

1*.0 

0.7086 

0.7*0.0 

0.7586 

0.7760 

0.7889 

0.7959 

0.8000 

i*.5 

0.721*2 

0.7760 

0.7937 

O.8139 

O.8181 

5.0 

0.7375 

0.7723 

0.7907 

0.8086 

0.8217 

O.829O 

0.8333 

5.5 

0.7^91 

O.803I* 

O.8213 

0.81*17 

«   0.8l*6l 

6.0 

0.7592 

0.7956 

O.81W* 

O.8323 

0.81*53 

O.8526 

0.8571 

6.5 

0.7682 

O.82U 

0.81*19 

0.8621 

0.8667 

7.0 

0.7763 

O.8138 

0.8327 

O.8501* 

O.8632 

O.870I* 

0.8750 

7.5 

0.7835 

0. 8U01* 

0.8579 

0.8777 

O.8823 

8.0 

0.7901 

O.8285 

0.81*73 

0.861*7 

0.8772 

0.881*2 

O.8889 

8.5 

0.7961 

O.8536 

O.8708 

O.8900 

0.891*7 

9.0 

0.8017 

0.81*07 

0.8593 

O.8763 

0.8881* 

0.8953 

0.9000 

9.5 

0.8068 

0.861*6 

O.881I* 

0.9002 

0.901*7 

10.0 

0.8115 

O.851O 

0.869!* 

0.8860 

O.8978 

0.901*5 

O.9091 

20 


TABLE  2 

PERIODIC 

-FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

cmin/cmax 

-  0.95 

(hA)« 

=  1.0 

1 

— — 

^r/^mln 

I 

NTU0 

1 

1.25 

1.5 

2 

3 

5 

00 

1 

O.U714 

0.^853 

0.1*901 

0.1*970 

0.5021 

0.501*7 

O.5063 

2 

0.6082 

0.6319 

0.61*5^ 

0.6592 

0.6693 

0.671*7 

0.6778 

3 

0.6760 

0.7061* 

0.7231* 

0.7^06 

0.7533 

O.7600 

0.761*0 

k 

0.7183 

0.7528 

0.7715 

0.7903 

0.8039 

O.8113 

O.8158 

5 

0.7^77 

0.781*9 

0.801*7 

0.821*0 

0.8379 

O.8U55 

O.8503 

6 

O.7698 

0.8090 

O.8292 

0.81*81* 

0.8621* 

0.8700 

0.8750 

7 

0.7872 

0.8277 

0.81*81 

0.8672 

0.8808 

0.8881* 

0.893^ 

8 

0.8012 

0.81*29 

O.8633 

0.8820 

0.8953 

0.9027 

0.9077 

9 

0.8129 

O.855I* 

0.8757 

0.891*0 

O.9069 

0.91^2 

0.9191 

10 

0.8229 

0.8661 

0.8861 

0.901*0 

0.9165 

0.9235 

0.9281* 

21 
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TABLE  3         PERIODIC-FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

Cmin/W  "     °-90 
(hA)«  =     1.0 


r "" cr/cmln 

1  1.25  1.5 


"1 

oo 


1  O.U763  0.4887  0.1+958  O.5030  O.5082  0.5110  0.5126 

2  0.6156  0.6405  O.65V7  0.6693  0.6800  O.6856  0.6889 

3  0.6846  0.7167  0.73^7  0.7530  0.7664  .     0.7735  0.7777 

4  0.7275  0.7642  0.7841  0.8040  O.8185  O.8263  O.83IO 

5  0.7573  0.7970  O.8181  O.8385  0.8534  0.8615  0.8664 

6  0.7795  0.8216  0.8432  6.8636  0.8784  O.8865  O.8916 

7  0.7970  0.8407  O.8625  O.8827  0.8972  O.9051  0.9102 

8  •  0.8112  O.8561  _  . 

O.8780  O.8978  0.9117  0.9195  0.9246 

9  O.8230  O.8689  O.8906  0.9099  0.9234  0.9309  0.9359 

10  O.8331  0.8797  0.9012  0.9199  0.9329  0.9401  0.9450 
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TABLE  4         PERIODIC-FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

(hA)«  =     1.0 


Cr/^n 

1.25  1.5 


oo 


1  0.4861  0.4991*  O.5071  O.5U9  O.5207  0.5254 

2  0.6299  0.6575  0.6750  0.6891  0.7011  0.7109 

3  O.7007  0.7364  0.7565  0.7769  0.7919  0.8045 
1+  0.7443  O.7853  0.8077  0.8299  0.8461  O.8597 

5  0.7743  O.819O  0.8427  O.8655  O.8818  O.8957 

6  0.7967  0.8440  O.8683  O.891O  O.9069  0.9206 

7  0.8140  0.8634  0.8879  0.9101  0.9253  0.9386 

8  0.8279  O.8788  0.9033  0.9249  0.9392  0.9518 

9  0.8394  0.8914  0.9157  O.9365  0.9500  O.9619 

10  0.8492  0.9021  0.9260  0.9460  0.9586  O.9696 
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TABLE  5         PERIODIC-FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

Cmin/Cmax  =     0.70 
(hA)«  «     1.0 


1 

cr/cmin 

NTU0 

1 

I.25 

1.5 

2 

3 

1 

0.1*959 

0.5104 

0.5187 

0.5271 

0.5333 

2 

0.6437 

0.6737 

0.6910 

0.7087 

0.7219 

3 

0.7155 

0.7548 

0.7770 

0.7995 

O.816O 

4 

0.7590 

0.8044 

0.8291 

0.8535 

0.8710 

5 

O.7887 

O.8382 

0.861*3 

O.8891 

O.9063 

6 

0.8104 

O.8629 

0.8895 

0.9139 

0.9303 

7 

O.8272 

O.8818 

O.9085 

0.9320 

0.9473 

8 

0.8405 

O.8967 

0.9232 

0.91*56 

0.9596 

00 

0.5384 
0.7326 
0.8295 
0.8855 
0.9207 


0.9439 
0.9598 
0.9709 

9  0.8514   0.9087   0.9348   0.9560   O.9687  0.9788 

10  0.8606    0.9188   0.9443    0.9642    0.9756  0.9845 
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TABLE  6         PERIODIC -FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

Cmin/Cfflax  =     0-50 
(hA)«  -     1.0 

I  ^r/^min   " — — 


1.25  1.5 


oo 


1  0.5156  0.5322  0.5^17  0.5515  0.5588  O.56U7 

2  0.6690  0.701*5  0.721*9  0.7^61  0.7618  0.77U6 
5  0.7^02  O.7869  0.8133  0.8399  0.8590  O.87M* 
1*  0.7817  O.8353  0.861*5  O.8928  0.9122  0.927^ 

5  O.8089  0.8670  0.8975  O.9256  0.9*06  0.9572 

6  0.8281*  0.8891*  0.9202  0.91*69  O.9629  0.971*5 

7  0.81*31  O.9062  O.9366  O.9615  0.9752  0.981*7 

8  0.851*6  0.9191  O.9U89  0.9717  O.9833  0.9908 

9  O.8639  0.9293  0.9582  0.9789  O.9886  0.99^ 

10  0.8716  0.9376  O.9655  0.981*2  0.9921  0.9966 
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TABLE  7 

PERIODIC 

-FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

Cmln/Cmax 

=  0.30 

(hA)« 

=  1.0 

1 

-  cr/cmln  - 

NTU0 

1 

1.25 

1.5 

2 

3 

1 

0.531*6 

0.5538 

0.561*8 

0.5761 

0.581*6 

2 

0.6902 

0.7315 

0.7556 

0.780U 

0.7988 

3 

0.7577 

0.811U 

0.81*22 

0.8729 

0.891*5 

1* 

0.7951 

0.8558 

0.8895 

0.921U 

0.9^20 

5 

O.8191 

O.8858 

0.9183 

0.9^90 

O.967O 

I 

00 

0.5915 
0.8136 
0.9110 
0.9566 
0.9787 

6  O.8360  O.9031  0.937^  O.9658  O.9808  O.9895 

7  0.81*87  0.9171  0.9506  0.9761*  O.9885  0.99^8 

8  O.8588  O.9278  O.9602  O.9835  0.9930  0.9971* 

9  0.8671  O.9363  0.9673  O.9881  0.9955  0.9987 

10  0.871*0  0.9*03  0.9729  0.9913  0.9971  0.9991* 
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TABLE  8    PERIODIC-FLOW  HEAT  EXCHANGER  EFFECTIVENESS 

CndnAW  ■  0.10 
(hA)'  -  1.0 

( cr/cmln 


1.25     1.5       2       3  oo 


1  0.5534  0.5752  0.5877  0.6008  0.6104  0.6186 

2  0.7070  0.7542  0.7820  O.8109  O.8319  0.8487 

3  0.7680  0.8279  O.863O  O.8981  0.9219  0.9391 

4  0.8009  0.8667  0.9042  0.9397  0.9615  0.9753 

5  0.8223  O.8908  O.9283  O.9618  0.9802  0.9900 

6  O.8378  0.9075  0.9440  0.9745  O.9892  0.9959 

7  0.8498  0.9200  0.9550  0.9824  0.9939  0.9983 

8  O.8595  O.9298  O.9632  0.9874  0.9964  0.9993 

9  0.8675  0.9376  0.9694  0.9908  0.9977  0.9997 

10  O.8743  0.9441  0.9743  0.9932  0.9984  0.9999 


27 


TABLE  9    PERIODIC -FLOW  HEAT  EXCHANGER  EFFECTIVENESS 
1 Cr/Cmln  


12        12        12 

NTU0  (hA)'  Cu&n/Cwax  "  1-0  Cmlri/Cwax  =  0-95  CmLn/Cmax   -  0.90 

3  1.00  0.6672  0.7280  0.6760  0.7406  0.6846  0.7550 

3  0.50  O.6676  0.7282  0.6752  0.7404  0.6826  0.7524 

3  0.25  0.6684  0.7284  0.6751  0.7404  O.6813  0.7521 

6  1.00  0.7592  O.8323  O.7698  0.8484  0.7795  O.8636 

6  O.50  0.7597  O.8323  O.7687  0.8481  0.7769  O.8628 

6  0.25  0.7608  O.832U  O.7685  0.8478  0.7751  0.8622 

9  1.00  0.8017  O.8763  0.8129  0.8940  0.8230  0.9099 

9  O.50  0.8021  O.8763  O.8115  O.8936  O.8198  0.9091 

9  0.25  0.8032  O.8763  0.8111  0.8933  O.8176  0.90811. 

Cmin/Cmax  -  °-Qo  cminAW  ■  0.70  0^/0^  =  0.50 

3  1.00  0.7007  0.7769  0.7155  0.7995  0.7402  0.8399 

3  0.50  O.6960  0.7752  0.7075  0.7964  0.7245  0.8331 

3  0.25  0.6919  0.7742  0.7001  0.79^3  0.7099  0.8270 

6  1.00  0.7967  0.8910  0.810U  0.9139  0.8284  0.9469 

6  O.50  0.7902  0.8888  0.8000  O.9103  0.8102  0.9398 

6  0.25  0.781*8  O.8872  0.7908  0.9074  0.7953  0.935^ 

9  1.00  0.8391*  0.9365  0.8514  0.9560  O.8639  O.9789 

9  0.50  O.8323  0.9344  0.8402  0.9529  0.8463  0.9735 

9  O.25  O.8263  0.9328  O.8306  O.9503  O.8327  0.9684 
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TABLE  9         CONTINUED 


( Cr/Cmin     ~ 1 

12  12 

NTU0     (hA)«           Cnjin/Cma^  =  0.30  Cmln/Cxnax  =  0.10 

3          1.00            0.7577        0.8729  0.7680        O.8981 

3          0.50            0.7338        0.8608  0.7373        O.8787 

3          0.25           0.7131        0.81*81*  0.7136        0.8579 

6          1.00           O.8360        O.9658  O.8379        0.971*5 

6          0.50            O.8130        0.95^6  O.8133        0.9597 

6          O.25            0.7957        0.9^39  0.7957        0.9^57 

9          1.00            0.8671        O.9881  0.8675        O.9908 

9          O.50            0.81*72        O.9808  0.81*72        O.9821 

9          O.25            0.8327        0.973^-  O.8326        0.9738 
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FIG.  13    AE/Eoo     vs  Cr/CjQin  in  the  range: 
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9.  Appendix  1  -  Derivation  of  Iterative  Equations 

Equations  [2]  ,  [3]  ,  [k]   and  [5J  concerning  the  element  on  the 
side  of  Cm<>Y,  Fig.  2(a),  from  the  body  of  the  thesis  are  repeated  here 
for  convenience. 


q  =  c,  (rrt  -TrO (■'/"*)        ■   ' 
q  --   ChA),  AT^?  ('M  'M)       .    . 


If  A^avg  from  equation  [5]  is  substituted  in  equation  (Uj  and  the 
value  for  Q  thus  obtained  is  substituted  in  equation  [2]   the  following 
will  be  obtained: 

This  may  be  solved  for  Txo,  giving 

T   -   f     <>*)»     1,  -r-    \    ^^"-^."l-r.      f/37 

In  a  similar  manner  Tro  may  be  solved  for  by  using  equation  /  3_/  *  giving 


'r"<3  ~ 


At  this  point  four  relations  can  be  derived  which  will  be  useful 
in  subsequent  manipulations  to  get  expressions  in  terms  of  the  dimen- 
sionless  parameters.  Prom  the  definition  of  NTU0  it  is  obvious  that: 


1*2 


An  expression  for  reciprocal  NTUX  may  be  derived  by  multiplying  and 
dividing  it  by  NTUn  as  expressed  in  equation  [15 j  t   thus, 

The  other  two ; expressions,  the  derivations  of  which  are  obvious,  are 
as  follows: 


<2v 


(fc/i)*  tfn/.  [/+ (m)'J 

and, 

Cr  CM)'  ^/L.* 


■       ■      p7j 


•  D*J 


To  obtain  TXQ  in  terms  of  only  the  inlet  temperatures  substitute 
equation  [l^J  in  equation  fl3j  and  solve  for  Tx0.  Starting  with  equa- 
tions [13]  and  [l^j  slightly  rearranged,  the  algebra  involved  follows. 


NTH* 


T*o    s    T~rL+     {to     + 


[fa-']*1  •  •  ■•  ■  -w 


-^  +  'JCr'i+^o)  +   1-   iML  +  )  \Tri    ...[H»-J 

ChA)* 


Substituting    [lW]   in   [13a]   , 


M<  '  1  J  J?>\ 


*3 


Each  of  the  terms  in  the  brackets  may  be  put  over  the  same  common  denom- 
inator which  may  then  be  multiplied  out,  giving 

Multiplying  through  by  (hA)'/2Nxcr  and  adding  and  subtracting  2Txi  on 
the  right  gives 


/V7K 


Now  tolving  for  Tx0  gives  a  form  convenient  for  iterative  calculations, 
thus, 

Tt0  ~-   Uc    -z,[>+   -7777;  +  —  j  tr,i  -Trl)   ,t..[i9] 

The  term  in  brackets  is  constant  for  a  given  problem  so  that 

T/0   -  -,/  -K,  (rti  -Tri)       .      ....     ,  ■  •  M 

where  Kx  may  be  put  directly  in  terms  of  the  dimensionless  parameters 

giving 

r         1  Jfo(hA)'  'I''  -  - 

To  obtain  Tro  in  terms  of  only  the  inlet  temperatures  substitute 
equatidn  [13]  in  equation  [l**-]  and  solve  for  TrQ.  Starting  with 
equations  [13]  and  \lkj   slightly  rearranged,  the  algebra  involved 
follows . 


1  a/7Wk    _/  L     ^    J 


1* 


/«'    "f    /<„      + 


AMt      A 
CM)*        I 


Tri 


Substituting   [iyo]  in    jito]   , 


(M)> 


*l 


J?A/iC 


■      -  KJ 


+   z  - 


K7U> 


Tti 


Multiplying  out  the  common  denominator  of  all  terms  in  brackets  gives 


-h 


4Jt  +  2  -zj  ra 


Multiplying  through  by  NTUx/2Nr  and  adding  and  subtracting  2Trl  on  the 
right  gives 


t 


(h/0*     /Vy-Cmo* 


+'    rr. 


^<c      ^  +,    rr£ 


t  2  (Tx<'  -7w) 


Now  solving  for  Tro  gives  a  form  convenient  for  iterative  calculations, 
thus, 


The  term  in  brackets  is  constant  for  a  given  problem  so  that 


ft6   *Trl    +  tfa  (r*<'  -  T"^3 


.  fclj 


^5 


■where  Kq   may  be  put  directly  in  terms  of  the  dimensionless  parameters 
giving 


K^'ZUt 


Qmm  Q x         V 


Cmax  Camo  Nt 


***  M  €Z™ 


M 


For  the  element  on  the  side  of  C^q  ,  Fig  3(h),  a  set  of  equations 
similar  to  equations  [s]  through  [5]  can  be  written.  They  are: 


Q  =  Cm;n  (J„,  -7n:)('/f«-)       ■    • 
0    --    C<   (Jri  -7rc){'/H<)  ■     ■ 

Q     -   (hA)n  A7I„a  ('A/,  /a/*)  .   . 

fat   -  Ch)(rr:  tTV.)  -  U^Oni  +  r«) 


.   .   .    .  [iQ 


Following  the  same  procedure  as  before  results  In  two  equations 
similar  to  equations  [13J  and  Jlk]   as  follows 


and 


t 


UA), 


*o  "/^c+Cm), 


^Cv,;„  +(hA)« 


2(V^Cr  -(KA)* 


•  f^J 


3AI„  C*  +  CM)n  _ 


n 


rt 


(T„i+Th0)  + 

These  may  be  solved  as  before  for  the  outlet  temperatures  in  terms  of 
only  the  inlet  temperatures  giving 


, .  .  I>] 


lm     -    Ini    -   l[_ 


.         ,  (-'»V\»VV  N'T 

1  f  —  Tu  * 


3  Mr 

V71U 


]"Vnc  -TTl)       '   "•    M 


k6 


and 


ro 


-    7Vi    +4>f    fc.ft    +     IMA    (^'"'0    ...f*J 


The  derivation  of  equations  [27]  and  [28  J  is  also  readily  appar- 
ent by  comparison  of  equation  [25]  with  [13]  and  [26J  with  [lU]  .  The 
only  difference  is  the  substitution  of  the  subscript  "min"  for  "max" 
and  "n"  for  "x"  and  then  equation  [27]  follows  from  [19J  and  equation 
f  28j  follows  from  [20J  . 

Equations  £27  J  and  [28J  may  be  written  as 


"^tio  ~  "nt  +  ^j  v '  fi   "    I ni  ) 


•     ■  f7j 


and 


/  VO   a   /y£   "~  f^f  C  '  *"'< '.    ~     I  ni  ) 


.  .  .  •  [7k] 


where 


f      1         2tfr         n 


-I 


and 


[tcj 


1   -£i-     / 

Equations  [13]  ,  [lU]  ,  (25  J  and  [26]  are  repeated  here  for  con- 


venience. On  the  side  of  C. 


max  ' 


T™  - 


Ch/i), 


«  UVrC^  +  CM)^ 


Uw  +  /ro)  +  / — TT7T     u 


r. 


ro 


i^fi-ii^lr^  ...f*7 


J^C<  +  IMJ* 


*7 


and  on  the  side  of  C^q  : 

T        ^  f  ("}& 1  , 

/no       [lMm;n  +LhA)J  &Yi+Tn)    + 


^  A/r  d;*  -  (mj> 


I  ro    ~ 


JtJfC„;„    t  (hA)n 


T*    ,  ,  .[*5] 


rt 


, .  N 


It  is  obvious  that  the  first  term  on  the  right  of  each  of  the  fore- 
going equations  is  a  positive  quantity.  It  can  next  be  observed  that 
the  coefficient  of  the  second  term  of  each  equation  must  also  be  positive 
for  if  it  were  not  the  outlet  temperature  in  question  would  grow  hotter 
as  its  inlet  temperature  became  colder,  or  vice  versa,  which  would  be 
contrary  to  the  second  law  of  thermodynamics. 

This  provides  a  convenient  set  of  conditions  to  be  met  before  the 
iterative  scheme  can  be  expected  to  converge.  By  using  the  relations  of 
equations  [15J  through  [l8J  these  conditions  may  be  stated  as  follows: 


HTUL0[l+    /W 


C/v\tf_X 


6  Z  t/r 


.      .      .19] 


UTU,[l  t    /(wj'j 


/ C  yyiLt 


Cr/n 


6  3Nt 


,  .  M 


A/TU 


[i-h   (Mj'j   ^  2^ 


.  [**] 


m0[n  (m)'J 


£     2*1  n 


,  -M 
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In  the  ranges  of  parameters  considered  and  by  using  an  equal  num- 
ber of  subdivisions  of  each  stream  it  can  be  seen  that  condition  ^9J  is 
the  most  severe  and  therefore  it  alone  can  be  used  as  the  condition  for 


convergence , 


k9 


10.  Appendix  2  -  Computer  Program  Details 

The  computer  used  for  this  work  was  the  National  Cash  Register  102A 
general  purpose  digital  computer  and  what  follows  assumes  the  reader  has 
some  familiarity  with  its  operation. 

The  basic  arrangement  of  the  computer  program  was  designed  to  accom- 
modate any  number  of  subdivisions  of  each  of  the  three  streams  up  to  and 
including  61*.  It  was  found  that,  on  the  average,  five  iterations  were 
required  to  meet  the  reversal  condition  and  limit  on  the  heat  balance 
error  with  a  Judicious  estimate  of  the  initial  matrix  temperature  dis- 
tribution. With  16  subdivisions  of  each  stream  the  time  required  for 
each  iteration  was  four  minutes.  For  different  degrees  of  subdivision 
the  time  required  for  each  iteration  would  vary,  as  would  be  expected, 
as  the  square  of  the  factor  by  which  the  subdivision  was  changed,  e.g. 
32  subdivisions  of  each  stream  would  require  16  minutes  for  each  iter- 
ation. It  was  determined  that  16  subdivisions  of  each  stream  could  be 
expected  to  give  accuracy  to  about  five  units  in  the  fourth  place  of  the 
effectiveness.  To  finish  in  a  reasonable  time  and  obtain  the  desired 
four  place  accuracy  it  was  decided  that  all  the  results  would  be  cal- 
culated using  16  subdivisions  of  each  stream  to  be  followed  by  suffic- 
ient check  calculations  with  32  subdivisions  to  permit  extrapolation  of 
the  results  to  infinite  subdivision. 

The  use  of  the  various  cells  of  the  memory  was  as  follows: 

Cells  Use 

0000-0^77  Program 

0500-0577  Initial  Tr's 

0600-0777         -  Program 

1000-1077  Current  parameters  and  K's,  and 

temporary  storage. 
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1100-1177  Outlet  temperatures  on  side  of  C^x 

1300-1377         Outlet  temperatures  on  side  of  Cn^n 
1500-1577  Initial-Intermediate-Final  Tr's 

The  composite  computer  program  flow  diagram  is  shovn  in  Fig.  Ik   and 
the  contents  of  cells  0000-0477  and  0600-0777  are  recorded  in  Table  10. 

Box  10  and  11  of  Fig.  lU  accomplish  the  iteration  process.  This  is 
done  by  the  repetitive  application  of  equations  [6J  to  each  element  on 
the  side  of  C^qx  in  box  10  and  equations  [7]  to  each  element  on  the  side 
of  Cjnin  in  box  11.  Since  the  internal  temperatures  are  of  no  concern 
in  this  problem  it  was  only  necessary  to  provide  permanent  storage  for 
the  fluid  outlet  temperatures  of  each  column  and  the  initial  and  final 
rotor  temperatures  of  each  row. 

Box  12  determines  if  the  reversal  condition  is  met  by  taking  the 
difference  between  the  initial  and  final  Tr  of  each  row.  A  difference 
greater  than  the  specified  limit  will  cause  the  iteration  cycle  to 
start  over.  Also  the  largest  difference  is  selected  so  that  if  the  re- 
versal condition  is  met  a  means  of  manual  inspection  is  provided,  if 
desired,  by  causing  this  difference  to  be  printed  out  upon  putting  switch 
2020  in  the  up  position  (box  17). 

If  the  reversal  condition  is  met  then  the  heat  balance  error  is 
computed  in  box  16  and  if  this  is  greater  than  the  specified  limit  the 
iteration  cycle  is  started  over.  The  equation  used  to  compute  the 
heat  balance  error  was  arrived  at  as  follows: 


n 
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where  Tno  and  Txo  are  the  average  fluid  outlet  temperatures.  Since 
TXi  was  taken  as  unity  and  Tn^  as  zero  then  the  heat  balance  error  re- 
duces to 


error 

V 


The  average  fluid  outlet  temperature  on  each  side  was  computed  by  sum- 
ming the  individual  outlet  temperatures  of  each  column  and  dividing  by 
the  number  of  subdivisions. 

Another  means  of  optional  manual  inspection  was  provided  in  box  19. 
If  desired,  the  heat  balance  error  would  be  converted  from  octal  to  decimal 
and  printed  out  by  putting  switch  20*^0  in  the  up  position. 

By  definition,  the  effectiveness  may  be  computed  by  either  of  the 
following  equations: 

c*(Thl  -  rhQ)  r 

t   * *"      .   .       ,   ..    (    ,  ,  [3laJ 

F  -  Cq   C^°  -  ^ ) •  OlloJ 

("ml*  {lUi    ~    Ui  ) 

where  T^q   and  Tco  are  the  average  fluid  outlet  temperatures  on  the  hot 
and  cold  sides,  respectively.  From  the  conditions  assumed  for  calcula- 
tion purposes  the  following  may  be  observed: 

Tcl  =  Tnl  =  0 

Thi  =  Txl  =  1 

cc  =  cmln  t   and, 

Tco  =  Tno 
and  therefore  the  expression  for  effectiveness  from  equation  plbj  re- 
duces to 

E  -  (Tno)avg       [32] 
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(Tno)aVg  will  have  already  been  computed  in  the  process  of  calculating 
the  heat  balance  error  and  therefore  in  box  21  it  is  only  necessary  to 
convert  this  to  decimal  and  print  it  out. 

The  final  Tr's  punched  on  cards  in  box  22  were  used  as  estimates 
for  the  initial  Tr's  in  subsequent  calculations.  This  was  usually  done 
by  using  the  final  Tr's  obtained  at  a  given  Cmin/Cmax  as  the  initial 
Tr's  for  the  next  lower  value  of  CWn/Cjaax,  at  corresponding  values  of 
the  other  parameters. 

The  remainder  of  the  program  consists  of  a  provision  for  interupting 
the  calculations  at  logical  stopping  places  and  a  system  of  inserting 
new  values  of  parameters  after  each  solution.  These  are  self-evident 
in  Fig.  14.  Parameters  are  entered  as  floating  point  numbers  and  the 
arithmetic  involved  in  calculating  the  K's  is  done  in  floating  point. 
The  K's,  scaled  down  by  a  factor  of  two,  are  then  converted  to  single 
precision  and  the  subsequent  iteration  is  accomplished  in  fixed  point, 
single  precision. 

The  basic  computer  program  which  is  recorded  in  Table  10  is  set 
up  to  use  16  subdivisions  of  each  stream.  To  modify  it  for  8  or  32 
subdivisions  the  following  cells  should  be  changed  as  indicated: 
Cell     8  Subdivisions     32  Subdivisions 

2105002100001*0 

7 
3505^0210015^0 
36200015^02001 
352000210011^0 
36153720002001 
36153720031537 


0307 

21050021000010 

0403 

5 

0406 

35051021001510 

0410 

36200015102001 

04ll 

352OOO2IOOIIIO 

0424 

36150720002001 

0425 

36150720031507 
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01*30 

35200021001310 

OU33 

2ff0003 

01*36 

30131001*332007 

01*37 

30111001*332007 

01*1*3 

35151021000510 

041*7 

36051015102000 

01*56 

01*300030000510 

0462 

05300030001510 

3520002100131*0 

2ff0005 
30131*001*332007 
30111*001*332007 
35151*02100051*0 
36051*0151*02000 
01*30003000051*0 
0530003000151*0 

In  its  present  form  the  program  by-passes  box  25  and  goes  directly 
to  box  27.  This  causes  all  calculations  to  be  made  with  (hA)'  =  1.  To 
get  results  with  (hA)'  =  l,  l/2,  and  l/l*  enter  the  following  changes: 
0007     21100601*020002 
O367     33100606760371 
0570     3^300021000373 

0371  361OO60UOOIOO6 

0372  31*300021000006 

In  its  present  form  the  program  produces  results  with  Cr/Cj^  = 
1,  1.25,  1.5,  2,  3,  and  5.  To  eliminate  Cy/C^n  =  5  enter  the  following 
changes : 

0375  3^067210020600 

0376  31*067310030600 

Space  is  provided  in  the  iteration  subroutines  for  entering  appro- 
priate commands  to  get  any  desired  temperatures  printed  out  before  they 
are  lost.  These  may  be  printed  in  decimal  simply  by  placing  in  cell 
2000  and  unconditional  transfer  to  the  decimal  conversion  and  print  sub- 
routine at  cell  0700. 


5^ 


BGK 

FUST 

COMMAND 

• 

1 

0000 

Set  current  Cmin/Cmax  ■  1 

2 

0001 

Set  current  Cr/Cmin  ■  1 

32 

0621 

Put  in  next  Cmin/Cmax 

/ 

3 

0002 

Set  current  NTUo  ■  1 

30 

0600 

Put  in  next  Cr/Cmin 

f 

k 

k 

0003 

Set  current  (hA)*  =  1 

28 

0603 

Put  in  next  NTUo 

* 

\, 

5 

0001* 

Print  current  parameters 

26 

0371 

Put  in  next  (hA)1 

6 

0010 

Compute  K's 

1 

7 

0126 

Test  switch  2100 
I 

No 
No 
No 
No 


10 


11 


12 


06U0 


0130 


013^ 


0200 


031U 


HL 


Down 


Fill  initial  Tr's  from  cards 


Initial  Tr's  double  stored 


Iteration  on  side  of  Gmax 


Iteration  on  side  of  Cain 


Reversal  condition  met? 


Compute 


15 


0306 


Print  final  Tr's  and  stop 
T 


UP 

L 


Down |ll+  j  030U~[  Test  switch  2010 


No 


->jl3 


Final  Tr's 


Initial 


Yes 
A. 


16    O2H0   I  Heat  balance  error  small  enough T 


No 


Yes 

^_^ 4 

I  17      j     0326        I  Test  switch  2020 

Down 

^_^ I  f 

I  19  0327        I  Test  switch  20^0 


Up 


»  18 


0331 


Print  largest  reversal 
temperature  difference 


T 


Up 


20 


0333 


Print  heat  balance  error 
in  decimal 


Down    , 

X     1 

21 

03^7 

Convert  E  to  decimal  and  print 

22 

0352 

Punch  final  Tr's  on  cards 

23 

036k 

Te8t  switch  2010 

1 

Up 


A 


0366 


atop 


Down 


Compute 


33 


31 


29 


27 


25 


0620 
0611 


0375 


0373 


0367 


Stop 

Cmin/Cmax  finished? 


Cr/Cmin  finished? 


NTUo  finished? 


(hA)'  finished? 
T 


Fig.  Ik       COMPUTER  PROGRAM  FLOW  DIAGRAM 
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TABLE  10 

0000  31040004011000 

0001  31040004011002 

0002  31040004011004 

0003  31040004011006 

0004  21100004020002 

0005  21100204020002 

0006  21100404020002 

0007  34300021000010 
0010  31040004012000 

0011  31100610072002 

0012  34300021000733 

0013  31200420052000 

0014  31100410052002 

0015  34300021000725 

0016  31200420052002 

0017  31040304012000 

0020  34300021000760 

0021  31200420051040 

0022  31200420052000 

0023  31100610072002 

0024  34300021000725 

0025  31200420051042 

0026  31100210032002 

0027  34300021000725 

0030  31200420051044 

0031  31040004012000 


COMPUTER  PROGRAM  COMMANDS 

0032  34300021000760 

0033  31200420051046 

0034  3H04210432000 

0035  34300021000725 

0036  312O0420O51050 

0037  3H00010012002 

0040  343OOO21000760 

0041  31200420051052 

0042  31100210032000 

0043  343O0O21O00725 

0044  31200420051054 

0045  31200420052002 

0046  31040004012000 

0047  34300021000760 

0050  31200420051056 

0051  31105210532000 

0052  31105610572002 

0053  34300021000733 

0054  31200420051060 

0055  31200420052002 

0056  31040004012000 

0057  34300021000733 

0060  31200420052002 

0061  34300021000760 

0062  31200420051010 

0063  3U05010512000 
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TABLE  10     CONTINUED 
0061*  31105^10552002 

0065  3^300021000733 

0066  31200U20052002 

0067  31040004012000 

0070  34300021000733 

0071  31200^20052002 

0072  3^300021000760 

0073  31200420051012 

0074  311040101*12000 

0075  31104610472002 

0076  34300021000733 

0077  31200420052002 
0100  31040004012000 

0101  34300021000733 

0102  312OO42OO52OO2 

0103  343O0O21OO0760 

0104  31200420051014 

0105  3H04410452000 

0106  31100210032002 

0107  34300021000733 

0110  31200420052002 

0111  31040004012000 

0112  34300021000733 

0113  31200420052002 

0114  34300021000760 

0115  31200420051016 


0116  30101110101011 

0117  30101310121013 

0120  30101510141015 

0121  30101710161017 

0122  33210010110010 

0123  33210010130010 

0124  33210010150010 

0125  33210010170010 

0126  17210030000640 

0127  3430OO210O013O 

0130  35040421000131 

0131  35052021001520 

0132  35013104050131 

0133  34040601310131 

0134  35041221000157 

0135  35041321000141 

0136  35041421000146 

0137  34300021000140 

0140  35040121002000 

0141  36200015202001 

0142  30200104002001 

0143  252OOHOII2OO2 

0144  25200110132003 

0145  36200020022000 

0146  35152020031520 

0147  34300021000154 
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TABLE  10  CONTINUED 

0150  00000000000000 

0202  35OI+252IOOO2II 

0151  00000000000000 

0203  35210021002000 

0152  00000000000000 

0201+  3611+7720002001 

0153  00000000000000 

0205  30200101+002001 

0151*  35011+101+07011+1 

0206  25200110152002 

0155  35011+601+05011+6 

0207  25200110172003 

0156  31*01+10011+1011+1 

0210  35200020022000 

0157  35200021001120 

0211  3611+77200311+77 

0160  31+300021000170 

0212  3^300021000216 

0161  00000000000000 

0213  00000000000000 

0162  00000000000000 

0211+  00000000000000 

OI63  00000000000000 

0215  00000000000000 

016U  00000000000000 

0216  360201+01+150201+ 

0165  00000000000000 

0217  36021101+050211 

0166  00000000000000 

0220  3^0201+01+260201+ 

0167  00000000000000 

0221  35200021001320 

0170  35015701+000157 

0222  3^300021000231 

0171  3^1101570135 

0223  00000000000000 

0172  31+300021000200 

0221+  00000000000000 

0173  00000000000000 

0 

0225  00000000000000 

0171+  00000000000000 

0226  00000000000000 

0175  00000000000000 

0227  00000000000000 

0176  00000000000000 

0230  00000000000000 

0177  00000000000000 

0231  35022101+000221 

0200  3501+2721000221 

0232  3to1<-3002210201 

0201  3501+21+21000201+ 

0233  31+300021000311+ 
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TABLE  10  CONTINUED 
0231*  00000000000000 

0235  00000000000000 

0236  00000000000000 

0237  00000000000000 
021*0  350U3^210002l*3 
021*1  350^3521000255 
021*2  35210021002000 
O2I4.3  30132001*332007 
021*1*  352OO02O072OO0 
021*5  35021*301*15021*3 
021*6  3l*Ol*3602l*302l*3 
021*7  30200001*00101*0 

0250  31210020002000 

0251  31100010012002 

0252  31*300021000725 

0253  312001*20052000 
0251*  3521O021002O03 

0255  30112001*332007 

0256  352OO32OO72OO3 

0257  35025501*150255 

0260  31*01*3702550255 

0261  3601*0120032003 

0262  31210020032002 

0263  35200201*002002 
0261*  31*300021000265 
0265  3^300021000760 


0266  3020052001*2005 

0267  3601*0120052005 

0270  30200501*002005 

0271  3^200501*1*00300 

0272  3^300021000326 

0273  00000000000000 
0271*  00000000000000 

0275  00000000000000 

0276  00000000000000 

0277  00000000000000 

0300  3501*1*221000300 

0301  35150021000500 

0302  35030101*050301 

0303  31*01*1*303010301 
0301*  17201030000306 
0305  31+300021000131* 

0306  2101*1*1*01*1*50001 

0307  21050021000020 

0310  22000000000000 

0311  31*300021000131* 

0312  00000000000000 

0313  00000000000000 

0311*  35210021001050 

0315  3501*1*621000316 

0316  36052015202000 

0317  31*200010500321* 
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TABLE  10  CONTINUED 

0320  35031604500316 

0321  340Wr03l603l6 

0322  34105004410300 

0323  34300021000240 

0324  35200021001050 

0325  34300021000320 

0326  17202030000331 

0327  17204030000333 

0330  343OO021OOO340 

0331  21105021000001 

0332  34300021000327 

0333  35200521002000 

0334  34300021000700 

0335  34300021000340 

0336  00000000000000 

0337  00000000000000 

0340  34300021000347 

0341  00000000000000 

0342  00000000000000 

0343  00000000000000 

0344  00000000000000 

0345  00000000000000 

0346  00000000000000 

0347  35104021002000 

0350  34300021000700 

0351  21044404450001 


0352  35045721000353 

0353  05300030001520 

0354  13046004550000 

0355  303OOOO4611077 

0356  13046OO455OOO4 

0357  30300004611077 

0360  35035304530353 

0361  34046203530353 

0362  34300021000363 

0363  34300021000364 

0364  17201030OOO366 

0365  34300021000367 

0366  22000000000000 

0367  34300021000373 

0370  00000000000000 

0371  00000000000000 

0372  00000000000000 

0373  34065210040603 

0374  34065310050603 

0375  34067410020600 

0376  34067510030600 

0377  34300021000611 

0400  00000000000000 

0401  00400000000000 

0402  02000000000000 

0403  00000000000006 
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TABLE  10     CONTINUED 

oi*oi*  35050021001500 

01*05  00000100000001 
Ol*06  35O5202IOOI52O 
01*07  00000000010000 

01*10  36200015202001 
01*11  35200021001120 
01*12  35200021001100 
01*13  36200015002001 
ol*ll*  35150020031500 

01*15  00000100000000 
Okl6  300500014.002000 
01*17  300520014.062000 
01*20  30152001*002000 

01*21  30150001*002000 
014-22  30110004002000 
014-23  3OII2OOI4.OO2OOO 
0l*2l*  36151720002001 
01*25  36151720031517 
01*26  3611*7720002001 
01*27  35200021001300 
01*30  35200021001320 
01*31  30130001*002000 
01*32  30132001*002000 
01*33  02000000000001* 
Ol*3l*  30130001*332007 
01*35  30110001*332007 


01*36  30132001*332007 
01*37  30112001*332007 
01*1*0  00000100000000 

01*1*1  00000001*000000 

01*1*2  35150021000500 
01*1*3  35152021000520 

01*1*4  0031*121231*1212 

01*1*5  10000000000000 
01*1*6  36050015002000 
01*1*7  36052015202000 
01*50  00000100010000 
01*51  OI+3OOO30000500 

01*52 15555555555556 

01*53  00000000000010 

01*51*  00000000000500 
01*55  U1 111.1111  nil 
01*56  01*300030000520 
01*57  053OOO30001500 
01*60  26666666666666 
01*61  00000000001000 

01*62  05300030001520 

01*63  00000200020000 
01*61*  310661*06651002 

01*65  31210006551000 

0U66  00000000000000 
01*67  00000000000000 
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TABLE  10  CONTINUED 
01+70  00000000000000 
01*71  00000000000000 
01+72  00000000000000 
01+73  00000000000000 
oW  00000000000000 
01+75  00000000000000 
01+76  00000000000000 
014-77  00000000000000 

0600  310661+06651002 

0601  35060001+630600 

0602  3^300021000002 

0603  3101+0001+012000 

0601+  311001+10052002 
0605  314.300021000733 

0606  312001+20051001+ 

0607  31040001+011006 

0610  3^300021000006 

0611  35OU6U210OO6OO 

0612  31100010012000 

0613  35210021002002 
0611*-  36210006632003 

0615  3^300021000733 

0616  3^200521000621 

0617  35046521000621 

0620  22000000000000 

0621  31210006551000 


0622  350621014O70621 

0623  3^300021000001 

06214-  00000000000000 

0625  00000000000000 

0626  00000000000000 

0627  00000000000000 

0630  00000000000000 

0631  00000000000000 

0632  00000000000000 

0633  00000000000000 
063*+  00000000000000 

0635  00000000000000 

0636  00000000000000 

0637  00000000000000 
061+0  350U5121000614-6 

061+1  05300030002100 

061+2  0601+5201+550000 

O6I+3  30300001+51+1077 

061+1+  0601+5201+550001+ 

061+5  30300001+5^1077 
061+6  01+300030000520 
061+7  35061+601+53061+6 

0650  3^01+56061+6061+1 

0651  31+300021000130 

0652  00000000000001+ 

0653  00500000000000 
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TABLE  10  CONTINUED 
0651*  00000000000001 

0655  0071*6311*6311*63 

0656  00711*6311*6311*6 

0657  0063114-6311*6315 

0660  0051+6311*6311*63 

0661  001*00000000000 

0662  002311*6311*6315 

0663  0006311*6311*632 

0661*  00000000000001 
0665  00500000000000 

0666  00000000000001 

0667  00600000000000 

0670  00000000000002 

0671  001*00000000000 

0672  00000000000002 

0673  00600000000000 
0671*  00000000000003 

0675  00500000000000 

0676  02000000000001 

0677  001*00000000000 

0700  27300007162001 

0701  350711*20010713 

0702  35200021000722 

0703  05300030000711* 

0701*  32200207002000 

0705  25200220072002 


0706  32200220012000 

0707  32210020012002 

0710  27200120052001 

0711  31*200121000705 

0712  21200007200001 

0713  31*300021000351 
0711*  31*300021000000 

0715  02000000000001* 

0716  02000000000030 

0717  00500000000000 

0720  01000000000000 

0721  0071*0000000000 

0722  001*63511721510 

0723  02000000000001 
0721*  00000000000001 

0725  30300007162001* 

0726  350711*2001*0732 

0727  35200020022001* 

0730  25200120032005 

0731  312001*20052001* 

0732  31*300021000253 

0733  30300007162001* 
0731*  350711*2001*071*5 

0735  3320002002071*6 

0736  36200020022001* 

0737  3020012001*2005 
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-   TABLE  10  CONTINUED 

071*0  35200321002006 

07^1  3520022100200** 
07^2  35200520062005 
07lv3  37200521000755 
OJkk   312001*20052001* 

07^5  3^300021000113 
071*6  362OO22OOO2OOU 
07U7  30200320042005 

0750  35200121002006 

0751  35200021002001* 

0752  31*50002100071*2 

0753  27200507232005 
0751*  30200507232005 

0755  2720050721*2005 

0756  352001*0721*2001* 

0757  31*30002100071*1* 

0760  3O3OOOO7I62OOI* 

0761  350711*2001*0766 

0762  36200020022001* 

0763  23200120032005 
0761*  37200521000767 

0765  312001*20052001* 

0766  31*300021000266 

0767  27200507232005 

0770  30200507232005 

0771  2720050721*2005 


0772  352001*0721*2004 

0773  31*300021000765 
0771*  00000000000000 

0775  00000000000000 

0776  00000000000000 

0777  00000000000000 
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